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(57) Abstract 



The invention relates to securing of information utilising optical imaging technologies and more specifically to phase encryption and 
decryption of images. An image is encrypted into a mask having a plurality of mask resolution elements (X m , Y m ) by encoding the image 
using e.g. a phase mask with an encoded phase value <f> (X m , Y m ) and an encoded amplitude value a (X m , Y m ), and by further encrypting 
the mask (using e.g. a spatial light modulator) by addition of an encrypting phase value <f> c (X m , Y m ) to the encoded phase value <f> (X m , Y m ) 
and by multiplication of an encrypting amplitude value ac (X m , Y m ) with the encoded phase value a (X m , Y m ). The method of decrypting 
comprises the steps of decrypting the mask by radiating electromagnetic radiation towards the mask and inserting into the path of the 
electromagnetic radiation a complex spatial electromagnetic radiation modulator comprising modulator resolution elements, the decrypting 
phase value <j>t (Xd, Yd) and the decrypting amplitude value ad (Xd, Yd) respectively, of a modulator resolution element (Xd, Yd) being 
substantially equal to -h£ c (X m , Y ra ) and ac" 1 (X m , Y m ). 
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AN OPTICAL ENCRYPTION AND DECRYPTION METHOD AND SYSTEM 

FIELD OF THE INVENTION 

5 The invention relates to securing of information utilising imaging technologies and more 
specifically to phase and amplitude encryption and decryption of images. 

BACKGROUND OF THE INVENTION 

10 It is well known to form an image by phase contrast imaging methods in which phase 
modulation of light is converted into intensity modulation. As opposed to intensity 
modulation, phase modulation does not involve loss of energy. 

In published patent application no. WO 96/34307, which is hereby incorporated by 
15 reference, a phase contrast imaging method is disclosed for calculating phasor values of a 
phase mask for synthesising a desired intensity pattern. 

In "Securing Information with Optical Technologies", Bahram Javidi, Physics Today, Vol. 
50, No. 3, March 1997, pp 27-32, a method and system for optically securing information is 

20 proposed. An image is encrypted in a 4f-lens optical configuration (i.e. comprising two 
Fourier transforming lenses) by Fourier transforming the image and an input phase mask 
with a first lens of the 4f-lens system. A Fourier plane phase mask having phasor values 
e ib(o,p) js p OS jtj onec | j n the Fourier plane of the first lens and an encrypted image is formed 
with a second lens that Fourier transforms images in the Fourier plane of the first lens. The 

25 encrypted image is decrypted in a similar 4Mens configuration in which another Fourier 
plane phase mask (the key) having phasor values e AbM) is positioned in the Fourier plane 
of the first lens. 

It is a disadvantage of the known cryptographic method that encrypting an image both in 
30 the object plane and the Fourier plane leads to generation of speckle patterns in the 
decrypted image thereby corrupting data having been encrypted. 

It is another disadvantage of the known cryptographic method that extremely accurate 
three-dimensional positioning of the phase mask in the Fourier plane is required for 
35 successful encryption and decryption. 
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It is yet another disadvantage of the known cryptographic method that both amplitude and 
phase have to be recorded in the encrypted mask. 

SUMMARY OF THE INVENTION 

5 

It is an object of the present invention to provide an apparatus of the above kind which 
apparatus is robust, compact, and simple to design and relatively cheap to manufacture. 

It is another object of the present invention to provide an improved method and apparatus 
10 without the above-mentioned disadvantages. 

It is still another object of the present invention to provide an improved method and 
apparatus for optically securing data utilising phase contrast imaging. 

15 According to a first aspect of the present invention, a method is provided of decryption of 
an encrypted image having a non-encrypted image intensity pattern l(x'y). 

According to a second aspect of the present invention, the image is encrypted into a mask 
having a plurality of mask resolution elements (x m ,y m ) by encoding the image into the mask 

20 with an encoded phase value <|>(x m ,y m ) and an encoded amplitude value a(x m ,y m ), and by 
encrypting the mask by addition of an encrypting phase value <j> c (x mj y m ) to the encoded 
phase value <(>(x m ,y m ) and by multiplication of an encrypting amplitude value a c (x m ,y m ) with 
the encoded amplitude value a(x mi y m ). Thus, each mask resolution element (x m ,y m ) 
modulates the phase and the amplitude of electromagnetic radiation incident upon it with 

25 the complex value a(x m ,y m )ac(x m ,y m )e^^ 

The method of decryption comprises the steps of decrypting the mask by radiating 
electromagnetic radiation towards the mask, and inserting into the path of the 
electromagnetic radiation a complex spatial electromagnetic radiation modulator 
30 comprising modulator resolution elements (x^), each modulator resolution element (Xd,y<j) 
modulating the phase and the amplitude of electromagnetic radiation incident upon it with a 
predetermined complex value ad(x dl y d )e l * d(Kd ' yd) 1 the decrypting phase value <j>d(Xd,y d ) and the 
decrypting amplitude value a d (Xd,y d ), respectively, of a modulator resolution element (x<j,y d ) 
being substantially equal to -<J>c(x mi ym) and a c " 1 (x m ,y m ), respectively, of a corresponding 
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mask resolution element (x m ,y m ), and imaging the mask and the electromagnetic radiation 
modulator onto the image having the image intensity pattern l(x\y'). 

Throughout the present description the combination of the mask and the electromagnetic 
5 radiation modulator is denoted an encoder. 

According to a third aspect of the present invention, a decryption system is provided for 
decrypting an encrypted image having a non-encrypted image intensity pattern l(x\y') that 
has been encoded into a mask having a plurality of mask resolution elements (x m? y m ) with 

10 an encoded phase value <|>(x m ,y m ) and an encoded amplitude value a(x m ,y m ), and encrypted 
by addition of an encrypting phase value <t> c (x m ,y m ) to the encoded phase value <Kx m ,y m ) and 
by multiplication of an encrypting amplitude value a c (x m ,y m ) with the encoded amplitude 
value a(x m ,y m ), each mask resolution element (x m ,y m ) modulating the phase and the 
amplitude of electromagnetic radiation incident upon it with the complex value 

15 a(x mi y m )a c (x m ,yje^ 

The system comprises a source of electromagnetic radiation for emission of 
electromagnetic radiation for illumination of the mask, a complex spatial electromagnetic 
radiation modulator that is positioned in the path of the electromagnetic radiation and 

20 comprises modulator resolution elements (x<j,y d ), each modulator resolution element (Xd,y d ) 
modulating the phase and the amplitude of electromagnetic radiation incident upon it with a 
predetermined complex value a d (x d ,y d )e i * d(xd yd) . The decrypting phase value <|) d (x d ,y d ) and 
the decrypting amplitude value a d (x d ,y d ), respectively, of a modulator resolution element 
(x^yd) are substantially equal to -<j> c (x m ,y m ) and a c * 1 (x m ,y m ), respectively, of a corresponding 

25 mask resolution element (x m ,y m ). 

Further, the system comprises an imaging system for imaging the mask and the 
electromagnetic radiation modulator onto the image having the image intensity pattern 

Kx'.y'). 

30 

The complex spatial electromagnetic radiation modulator may be positioned anywhere in 
the path of the electromagnetic radiation. When the complex modulator is positioned 
adjacent to the mask, the encrypted mask is decrypted when the decrypting phase value 
<MXd,yd) and the decrypting amplitude value a d (x d ,y d ), respectively, of a modulator 
35 resolution element (Xd,y d ) being substantially equal to -<j> c (Xm,ym) and a c ~ 1 (x m ,y m ), 
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respectively, of a corresponding optically aligned mask resolution element (x m ,y m ).The 
mask and the decrypting complex modulator are optically aligned when corresponding 
mask resolution elements (x m> y m ) and complex modulator resolution elements (Xd t y d ) are 
imaged onto the same desired resolution element (x\y') of the image having the image 
5 intensity pattern l(x\ y 1 ). 

For example, the complex modulator may be physically positioned remotely from the mask 
and imaged onto the position of the mask whereby the complex modulator is virtually 
positioned adjacent to the mask. Thus, the complex modulator may be positioned virtually 
10 or physically in another position than adjacent to the mask still performing the decrypting 
function when the amplitude and phase values of the complex modulator are transformed 
accordingly, e.g. by a Fresnel propagator. 

As already mentioned, the optically aligned combination of the mask and the decrypting 
15 complex modulator is denoted the encoder. The virtual mask resolution elements of the 
encoder are denoted (x,y) and thus for optically aligned resolution elements, each of 
resolution elements (x,y) and (x m ,y m ) and (Xd,y d ) is imaged onto the same image resolution 
element (x'.y 1 ). For example when the mask and the complex modulator are positioned 
adjacent each other: 

20 

(x,y) = (x m ,y m ) = (x d ,y d ). 

It is an important advantage of the present invention that the encrypted image may 
selectively require recording in the mask of phase values or amplitude values or a 
25 combination thereof. 

It is another advantage of the present invention that decryption is performed in a plane 
adjacent to the mask or an equivalent plane whereby generation of speckles in the 
decrypted image is suppressed. 

30 

It is yet another advantage of the present invention that there is no requirement of 
positioning a radiation modulator in the Fourier plane whereby an accurate three- 
dimensional positioning requirement is avoided. 
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It is still another advantage of the present invention that only a single key for encryption is 
required. 

It is preferred to encode the non-encrypted image into a phase mask and encrypt the 
5 phase mask by adding encrypting phase values to the encoded phase values. An 
encrypted phase mask is extremely difficult - if not impossible - to replicate by 
counterfeiters. Further, phase masks may be readily produced while masks requiring 
recording of both amplitude and phase are extremely complicated to produce, typically 
requiring production of two masks to be accurately superpositioned. The same applies to 
10 the decrypting complex modulator. 

The imaging may be performed with a common path interferometer, such as a phase 
contrast imaging system, a dark field imaging system, a field absorption imaging system, a 
point diffraction imaging system, a Smartt interferometer, a Schlieren interferometer, etc, 
15 and any combination hereof. 

Thus, the imaging for reconstruction of the intensity pattern l(x\y') may comprise: 

Spatial modulation of electromagnetic radiation with the encoder for modulation of the 
20 phase of the incident electromagnetic radiation by phasor values of individual resolution 
elements of the encoder, each phasor value e l,Kx y) being determined in such a way that 

1) the values of the Fourier transformed phasors attains predetermined values for 
predetermined spatial frequencies, and 
25 2) the phasor value of a specific resolution element of the encoder corresponds to a 
distinct intensity level of the image of the resolution element in the intensity pattern, 

Spatial phase filtering of electromagnetic radiation with a spatial phase filter for phase 
shifting of a part of the electromagnetic radiation, and 

30 

Generating the intensity pattern l(x\y') by interference in the image plane of an imaging 
system between the part of the electromagnetic radiation that has been phase shifted by 
the phase filter and the remaining part of the electromagnetic radiation. 
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Preferably, the method of decryption further comprises the steps of, and the decryption 
system further comprises means for Fourier or Fresnel transforming electromagnetic 
radiation modulated by the mask and the complex spatial electromagnetic radiation 
modulator, Filtering the Fourier or Fresnel transformed electromagnetic radiation by, in a 
5 region of spatial frequencies comprising DC in the Fourier or Fresnel plane, phase shifting 
with a predetermined phase shift value 6 the modulated electromagnetic radiation in 
relation to the remaining part of the electromagnetic radiation, and multiplying the 
amplitude of the modulated electromagnetic radiation with a constant B, and in a region of 
remaining spatial frequencies in the Fourier or Fresnel plane, multiplying the amplitude of 
10 the modulated electromagnetic radiation with a constant A, forming the intensity pattern by 
Fourier or Fresnel transforming, respectively, the phase shifted Fourier or Fresnel 
transformed modulated electromagnetic radiation, whereby each resolution element (x m ,y m ) 
of the mask is imaged on a corresponding resolution element (x\y') of the image, the 
filtering parameters A, B, e substantially fulfilling the equation 

15 

l(x\y') = A 2 |a(x , ,y , )e l * (x ' y ' ) + a (BA" 1 e l9 -1)| 2 

a being the average of the complex phasors a(x,y)e i<Kxy) of the encoder. 

20 It should be noted that, in each resolution element of the encoder, one of two phasor 
values which represent a particular grey level of the intensity pattern Kx'Y) may be 
selected. 

It is further preferred that the filtering parameters A and B substantially fulfil that A=1 and 
25 B=1. 

Preferably, the absorption of the mask is substantially uniform. 

It is also preferred that the phase shift value 6 substantially fulfils the equation 

30 



Ci— ! 



2|sin|| 
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and in this case 



5 /(x 1 ) y) = 2^1 + sin((t > --^(x 1 ,y) + |) 

For selected phase shift values <|>, <|>- being the phase of a . 

In a preferred embodiment of the present invention, the phase shift 6 is substantially equal 

10 to 7i. According the previous equation 9 = n leads to |a| = Yt, and then the phasor values 
e *(x.y) 0 f tne enc0C | er may De calculated in accordance with 

l(x',y') = 2(1-cos(<Kx\y'))) 

15 J ^sm^x^dxdy = 0 . 

encoder 

Although, the present methods of encryption and decryption are related to encoding in two 
spatial dimensions (planar encoding), the principles of the methods may be utilised for 
20 encoding in one to three spatial dimensions and/or in the temporal dimension. 

The electromagnetic radiation may be of any frequency range of the electromagnetic 
spectrum, i.e. the gamma frequency range, the ultraviolet range, the visible range, the 
infrared range, the far infrared range, the X-ray range, the microwave range, the HF (high 
25 frequency) range, etc. The present invention is also applicable to particle radiation, such as 
electron radiation, neutron radiation, etc. 

Preferably, the electromagnetic radiation is monochromatic or quasi-monochromatic so that 
the energy of the electromagnetic radiation is concentrated in a narrow frequency 
30 bandwidth. Since the intensity pattern is reconstructed by interference of two 

electromagnetic waves emitted from a common source of electromagnetic radiation, the 
phases of which have been changed differently, it is required that the frequency range of 
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the emitted electromagnetic radiation is sufficiently narrow to ensure that the two waves of 
electromagnetic radiation are coherent so that their superposition generates the desired 
intensity pattern. If the frequency range is too broad, the two waves will be incoherent and 
the phase information will be lost as superposition of non-coherent waves results in a 
5 summation of the intensities of the two waves. It is required that the difference between 
individual delays of electromagnetic radiation to be superpositioned is less than the 
wavelength of the radiation. This is a relaxed requirement that allows the electromagnetic 
radiation to be relatively broad-banded. For example in the visible range a Xe-lamp or a 
Hg-lamp can be used as a light source in a system according to the present invention with 
10 the advantage compared to a laser light source that speckle noise is reduced. The 
requirements of the spatial coherence of the electromagnetic radiation depend upon the 
space bandwidth product of the corresponding system and how close the required system 
performance is to the theoretically obtainable performance of the system. 

15 Preferably, the electromagnetic radiation is generated by a coherent source of 

electromagnetic radiation, such as a laser, a maser, a phaselocked laser diode array, etc. 
However a high pressure arc lamp, such as a Hg lamp, a Xe lamp, etc, may also be used 
and even an incandescent lamp may be used as a source of electromagnetic radiation in a 
low performance system. 

20 

The encoder changes the phase of an electromagnetic wave incident upon it. The encoder 
may transmit or reflect the incident electromagnetic wave. The encoder is divided into a 
number of resolution elements, each of which modulates the incident electromagnetic wave 
by changing its phase by a specific predetermined value. The predetermined values are 

25 assigned to each resolution element in different ways depending upon the technology 
applied in the component. For example in spatial light modulators, each resolution element 
may be addressed either optically or electrically. The electrical addressing technique 
resembles the addressing technique of solid-state memories in that each resolution 
element can be addressed through electronic circuitry to receive a control signal 

30 corresponding to the phase change to be generated by the addressed resolution element. 
The optical addressing technique addresses each resolution element by pointing a light 
beam on it, the intensity of the light beam corresponding to the phase change to be 
generated by the resolution element illuminated by the light beam. 
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Spatial phase modulation may be realised utilising fixed phase masks, devices comprising 
liquid crystals and being based on liquid crystal display technology, dynamic mirror 
devices, digital micro-mirror arrays, deformable mirror devices, membrane spatial light 
modulators, laser diode arrays (integrated light source and phase modulator), smart pixel 
5 arrays, etc. 

A spatial phase filter is typically a fixed phase mask, such as an optically flat glass plate 
coated with a dielectric layer at specific positions of the glass plate. However, the spatial 
phase modulators mentioned in the previous section may also be used for spatial phase 
10 filters. 

The imaging system maps the phase modulating resolution elements of the encoder on the 
target surface of the reconstructed intensity pattern. It may comprise a 4f-lens configuration 
(two Fourier transforming lenses utilising transmission of light or one Fourier transforming 
15 lens utilising reflection of light) or a single imaging lens. However, any optical imaging 
system providing a filtering plane for the spatial phase filter may be applied in a phase 
contrast imaging system. 

In the method of decryption according to the present invention, the reconstructed intensity 
20 pattern is generated by superposition of two electromagnetic waves in the image plane of 

the imaging system. The encoder changes the phase values of an electromagnetic wave 

incident upon it and the imaging system directs the electromagnetic wave with changed 

phases reflected from or transmitted through the encoder towards the spatial phase filter. 

The phase filter phase shifts a part of the electromagnetic radiation and the imaging 
25 system is adapted to superimpose in the image plane the phase shifted part of the 

electromagnetic radiation with the part of the electromagnetic radiation that is not phase 

shifted by the spatial phase filter. 

According to a preferred embodiment of the invention, the encoder is positioned at the front 
30 focal plane of a lens while the spatial phase filter is positioned in the back focal plane of the 
lens, whereby a first electromagnetic field at the encoder is Fourier transformed by the lens 
into a second electromagnetic field at the phase filter. Thus, specific spatial frequencies of 
the first electromagnetic field will be transmitted through the spatial phase filter at specific 
positions of the phase filter. For instance, the energy of the electromagnetic radiation at 
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zero frequency (DC) is transmitted through the phase filter at the intersecting point of the 
Fourier plane and the optical axis of the lens also denoted the zero-order diffraction region. 

It is presently preferred that the spatial phase filter is adapted to phase shift the DC-part of 
5 the electromagnetic radiation and to leave the remaining part of the electromagnetic 
radiation unchanged or, alternatively, to leave the DC-part of the electromagnetic radiation 
unchanged and to phase shift the remaining part of the electromagnetic radiation. The last 
alternative is preferred when the energy level of the DC-part of the electromagnetic 
radiation is so high that the phase shifting part of the phase filter will be destroyed by it. For 
10 example in laser cutting, the DC level of the laser beam can be so high that a phase 
shifting dot positioned at the intersecting point of the DC part of the laser beam at the 
phase filter would evaporate. It is also possible to block the electromagnetic radiation (no 
transmittance) in the zero-order diffraction region, however, the DC energy of the radiation 
is then lost. 

15 

Below, an expression that is valid for phase modulation of the intensity of the reconstructed 
intensity pattern as a function of the phasor values e l<Mx,y) of the encoder, when the DC-part 
of the electromagnetic radiation is phase shifted, is deduced. 

20 Electromagnetic radiation incident on the encoder can be described by a function A(x,y), 
where A(x,y) is a complex number (amplitude and phase) of the incident field on the 
resolution element (x,y) of the encoder. At the resolution element (x,y), the encoder 
modulates the phase of the incident radiation with a value c|>(x,y) so that the field after 
reflection by or transmission through the encoder may be described by the function A(x,y) * 

25 e** (x y) , e'* (x,y) being the phasor value of the resolution element (x,y) of the encoder. As A(x.y) 
preferably is a constant value over the entire surface of the encoder, the term is left out of 
the following equations for simplicity. 

The expression of the electromagnetic radiation incident on the spatial phase filter may 
30 now be separated into an AC-term and a DC-term. If the DC-term of the field is denoted a , 
the AC-term of the field is given by the term e'" ix ' y) - a . Since the spatial phase filter 
changes the phase of the DC-part of the electromagnetic radiation by 8, the intensity of the 
reconstructed intensity pattern at the image plane of the imaging system is given by: 



35 



l(x',y')= |e^+a(e ie -1)| 2 
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wherein (x\y') is the co-ordinates of the image of the resolution element (x,y) of the 
encoder formed by the imaging system in the image plane. 

5 It should be noted that the second term of the equation is a complex number that adds to 
the phasors e* (x,y) of the encoder and may be interpreted as a contrast control parameter 
for the reconstructed intensity pattern l(x',y'). 

According to a preferred embodiment of the invention, the average value of the phasors is 
10 adjusted in order to control the range of intensity levels. 

Instead of phase shifting the DC-part of the electromagnetic radiation, it is also possible to 
synthesise a prescribed intensity pattern by phase shifting other parts of the 
electromagnetic radiation by adapting the spatial phase filter to phase shift electromagnetic 

15 radiation incident upon one or more arbitrary regions of the phase filter and leaving the 
phase of the remaining part of the electromagnetic radiation unchanged and then 
superimposing the two parts of the electromagnetic radiation. The corresponding 
mathematics and the corresponding design procedures for the encoder and spatial phase 
filter will of course be more complicated than for the method described in the previous 

20 section. 

A simple example of phase shifting a part of the electromagnetic radiation of a spatial 
frequency different from the zero frequency is provided by moving the DC-part of the 
electromagnetic radiation to another spatial frequency in the Fourier plane (identical to the 
25 plane of the spatial phase filter) utilising an optical component with an appropriate carrier 
frequency (i.e. a grating or a prism) or, preferably, encoding the function of a grating or a 
prism into the encoder, and adapting the spatial phase filter to change the phase of the 
electromagnetic radiation at this spatial frequency and to leave the phase of the remaining 
part of the electromagnetic radiation unchanged. 

30 

According to another preferred embodiment of the invention, the encoder is not positioned 
in the back focal plane of the lens but in the Fresnel region of the lens instead. In this case, 
the electromagnetic field at the phase filter will be given by a Fresnel transformation of the 
electromagnetic field at the encoder. This further complicates the mathematics and the 
35 design procedures, for example the term in equation (7) below has to be substituted by the 
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value of the Fresnel transformation at the point(s) of phase changes of the phase filter. 
However, the Fresnel transformation may be calculated from a Fourier transformation by 
multiplication of the phasor values of the encoder by a quadratic phase factor followed by a 
Fourier transformation. 

5 

It is an important aspect of the present invention that each intensity level of the 
reconstructed intensity pattern for each resolution element may be generated by at least 
two different phasor values of a resolution element of the encoder. 

10 For example, when the spatial phase filter phase shifts the DC-part of the electromagnetic 

radiation, it will be shown later that, advantageously, the average a of the phasors of the 
resolution elements of the encoder should be equal to V2 and the value of the phase shift 9 
should be equal to n. In this case, the intensity of the reconstructed image pattern at the 
image (x\y ! ) of the resolution element (x,y) will be given by: 

15 

i(x\ y) = 2(i-cos<t)(x j y)) 

It is seen that complex conjugate phasors (values of <j> of opposite sign) result in identical 
intensity levels l(x\y'). It can be shown that for any value of the modulus of the average of 
20 the phasors |a| two phasors exist that will generate identical intensity levels of the 
reconstructed intensity pattern. 

Further, if the spatial phase filter phase shifts parts of the electromagnetic radiation 
different from the DC-part, the phasor value that generates a specific intensity level will 

25 depend on the position of the resolution element in question, i.e. the phasor value and the 
position of the resolution element with that phasor value together define the intensity level 
at the image of the resolution element in the reconstructed intensity pattern. Still, it is true 
that for each resolution element of the encoder, each intensity level of the reconstructed 
intensity pattern may be represented by one of two different phasors of complementary 

30 phase values. 

This freedom of being able to select, for each intensity level to be generated and for each 
resolution element of the encoder, one of two phasors is used to control the phase of the 
Fourier transform of the phasors at specific spatial frequencies by selection of phasors with 
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appropriate phase values to ensure two intervals of biunique functional dependence 
between phasor values and corresponding intensity values. 

This freedom of choice of phasors may be utilised to select phasors of neighbouring 
5 resolution elements of the encoder with a maximum difference between them, thereby 
generating an electromagnetic radiation emitted from the encoder with a maximum content 
of high spatial frequencies which will generate a good separation of the DC part of the 
electromagnetic radiation from its AC part. However, any other strategy of selecting 
between two possible phasor values of each resolution element may be chosen to 
10 generate a desired spatial frequency content of the electromagnetic radiation. 

Preferably, the phase of the Fourier transform of the phasors at specific spatial frequencies 
is adjusted in order to control whether the relation between each phasor and the 
corresponding intensity level is a monotonic increasing or a monotonic decreasing function. 

15 

Below, a set of different steps that may be incorporated - either alone or in combination - 
into the method of encryption is described. They are provided according to the present 
invention for adjustment of the modulus of the Fourier transform of the phasors at specific 
spatial frequencies to attain a prescribed value. 
20 According to one of the steps, the individual phasors of the resolution elements of the 
encoder are adjusted by a constant value until the desired value of the modulus of the 
Fourier transform of the phasors at specific spatial frequencies is attained while 
maintaining prescribed relative intensity levels between intensities of resolution elements of 
the intensity pattern, i.e. iteratively. 

25 

According to another step, the individual phasors of the resolution elements of the encoder 
are adjusted utilising histogram techniques known from image processing. A histogram is a 
bar chart showing the number of resolution elements of the reconstructed intensity pattern 
with a specific intensity value as a function of the intensity value. Any histogram technique, 
30 such as histogram equalisation, adapting the histogram to a predetermined distribution, 
etc., may be used iteratively until the modulus of the Fourier transform of the phasors at 
specific spatial frequencies attain the prescribed value. 
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According to yet another step, the phasor pattern of the encoder is spatially scaled in order 
to adjust the modulus of the Fourier transform of the phasors at specific spatial 
frequencies. 

5 According to still another step, the modulus of the Fourier transform of the phasors at 
specific spatial frequencies is adjusted utilising half tone coding techniques, such as raster 
techniques, area ratio modulation, spot diameter modulation, etc. 

Each complex phasor a(x,y)e l,Kx,y) of the encoder may be selected from a set of two 
10 determined phasors with complementary complex phasor values a(x,y)e i * 1(x,y) and 

a(x,y)e l,|,2(x ' y) in such a way that a specific spatial frequency distribution of the intensity of the 
electromagnetic radiation in the Fourier or Fresnel plane is attained. 

The phase <|>(x,y) of complex phasors a(x,y)e i<Kx,y) of adjacent resolution elements may 
15 alternate between the two possible complementary complex phasor values a(x,y)e i,{>1(x,y) and 

a(xy)e i*2(x.y) 

The complex phasors a(x,y)e i<|,1(xy) and a(x,y)e i<l>2(x ' y) may be complex conjugated. 

20 It is seen from the description above that the intensity levels may differ from one 
reconstructed intensity pattern to the next as a consequence of the adjustments of the 
modulus of the Fourier transform of the phasors at specific spatial frequencies. Thus, it is 
preferred to control the power of the radiation source in dependence of the intensity range 
of the intensity pattern so that a sequence of different intensity patterns shows uniform 

25 intensity levels. 

According to a preferred embodiment of the invention, the shape of the phase filter is 
adapted to match the spatial frequency content of the phasors of the encoder, e.g. to 
optimise the desired separation of the part of the electromagnetic radiation to be phase 
30 filtered from the remaining part of the electromagnetic radiation. 

It is within the scope of the present invention that the imaging system further comprises 
zooming means for variable scaling of the reconstructed intensity pattern. The zooming of 
the imaging system may be dynamically controllable, e.g. in response to the scaling of the 
35 pattern of phasor values of the encoder. 
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According to the present invention, the power of the radiation source may be controllable in 
response to the spatial scaling of the pattern in the encoder and/or the zooming of the 
focusing system. 

5 

In order to provide a compact and integrated system according to the present invention, the 
optical function of a Fourier-transforming lens is encoded into the phasors of the encoder. 
The Fourier transforming lens may be refractively or diffractively encoded. 

10 Similarly, the optical function of an output lens may be encoded into the phase filter either 
refractively or diffractively. 

Further, compensation may be encoded into the phasor values of the encoder so that part 
of the electromagnetic radiation modulated by the encoder has a substantially flat intensity 
15 profile in the image plane. Without this compensation, part of the electromagnetic radiation 
modulated by the encoder will have a flat profile with perturbations resulting from the phase 
filtering superpositioned upon it This may cause "ringings" (oscillations) at the edges of the 
reconstructed intensity pattern. 

20 According to another preferred embodiment of the invention, the source of electromagnetic 
radiation comprises one or more light sources of different wavelengths corresponding to 
three different colours, such as red, green and blue, for generation of intensity patterns of 
arbitrary colours. Further, several independent systems each one illuminated by its own 
wavelength can be combined into a single multi-wavelength system. 

25 

According to yet another preferred embodiment of the present invention, the phasor values 
of the resolution elements of the mask and the electromagnetic radiation modulator are 
binary, i.e. each phasor value is selected from a set constituted by two values, for example 
1 and -1 (corresponding to phase values 0 and n, respectively). It is further preferred that 
30 the random distribution of the decrypting binary phasor values is uniform. 

It is also preferred that the binary values are selected so that encryption and decryption 
correspond to a binary XOR logical operation since the XOR operation provides maximum 
ambiguity as to which inputs generate a specific output. 
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It is even more preferred that there is no spatial correlation between resolution elements in 
the encoded, non-encrypted phase pattern. 

The reconstructed intensity pattern l(x\y') may be detected by a camera, such as a CCD 
5 camera, and the recorded image may be transmitted to a computer for recognition of the 
image l(x',y*) and thus, authentication of the encrypted mask. 

BRIEF DESCRIPTION OF THE DRAWINGS 

10 Fig. 1 shows a 4f common path interferometer, 

Fig. 2 shows a 2f common path interferometer, 

Fig. 3 shows a 1f common path interferometer, 

15 

Fig. 4 shows (A) off-axis read-out of reflective SLM and (B) on-axis read-out of reflective 
SLM, 

Fig. 5 shows schematically an example of a prescribed intensity pattern in 1 D, 

20 

Fig. 6 shows schematically the resulting phase encoding corresponding to Fig. 5, 

Fig. 7 illustrates a binary phase image encryption method according to the present 
invention, 

25 

Fig. 8 illustrates a binary phase image decryption method according to the present 
invention, and 

Fig. 9 illustrates alignment of an encrypted phase mask and a decrypting phase mask. 

30 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

Fig. 1 shows a 4f common path interferometer 1 . A laser 2 emits a light beam which is 
expanded by a beam expander 3 into a plane light wave of uniform intensity and directs it 
35 towards an encoder 4, i.e. a combination of an encrypted mask 4a and a decrypting 
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complex spatial electromagnetic radiation modulator 4b. The light beam is transmitted 
through the encoder 4 and a Fourier transforming lens 5. The encoder 4 is positioned in the 
front focal plane of the lens 5 and a spatial filter 6 is positioned in the back focal plane of 
the lens 5 that is also the front focal plane of a lens 7. The Fourier transforming lenses 5, 7 
5 need not have identical focal lengths. Different focal lengths lead to a magnification ratio 
different from one. The filter 6 phase shifts by 9 and optionally attenuates (by a factor B) 
the zero order diffraction part 8 of the light modulated by the encoder 4. Optionally, the 
remaining diffraction part of the light modulated by the encoder 4 may be attenuated by a 
factor A. The reconstructed intensity pattern is generated in the back focal plane 9 of the 
10 lens 7 and a dynamic focusing system 10 images the reconstructed intensity pattern l(x\y') 
onto a focusing plane 1 1 . The reconstructed intensity pattern l(x\y') may be detected by a 
camera, such as a CCD camera, and the recorded image may be transmitted to a 
computer 12 for recognition of the image l(x\y') and thus, authentication of the encrypted 
mask 4a. 

15 

The computer 12 controls the optical system. The computer 12 comprises interface means 
for addressing each of the resolution elements of the complex spatial electromagnetic 
radiation modulator 4b and transmitting a decrypting value to the addressed resolution 
element. Alternatively, the complex spatial electromagnetic radiation modulator 4b may be 
20 a fixed mask i.e. a mask manufactured with fixed modulating values of the resolution 
elements. 

Optionally, the phase shift 0 and attenuation factors (A, B) of the filter 6 is adjustable and 
controllable by optional phase control means of the computer 12 which may be further 
25 adapted to adjust the phase shift, e.g. utilising equation 18. 

Fig. 2 shows a 2f common path interferometer 20. A laser 21 emits a light beam which is 
expanded by a beam expander 22 into a plane light wave of uniform intensity and directs it 
towards an encoder 23, i.e. a combination of an encrypted mask 23a and a complex spatial 

30 electromagnetic radiation modulator 23b and a polarisation beam splitter 24 and a quarter- 
wave plate 25. The polarisation beam splitter 24 and the quarter-wave plate 25 allows 
beam-splitting of light of a specific linear polarisation without the power loss associated 
with conventional beam-splitters due to splitting of the beam in both directions of 
transmission through the beam-splitter. After transmission through the polarisation beam 

35 splitter 24 and the quarter-wave plate 25, the light beam is transmitted through a Fourier 
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transforming lens 26 and is reflected from a spatial filter 27. The encoder 23 is positioned 
in the front focal plane of the lens 26 and the spatial filter 27 is positioned in the back focal 
plane of the lens 26. The filter 27 phase shifts by 0 and optionally attenuates (by a factor B) 
the zero order diffraction part 28 of the light modulated by the encoder 23. Optionally, the 
5 remaining diffraction part of the light modulated by the encoder 23 may be attenuated by a 
factor A. The reconstructed intensity pattern l(x\y') is generated in the back focal plane 29 
of the lens 26 and a dynamic focusing system 30 images the reconstructed intensity 
pattern l(x\y') onto a focusing plane 31. As described for the system shown in Fig. 1, the 
image l(x',y') may be detected with a camera and transmitted to a computer 32 for 
10 processing and authentication and the system 20 may also be controlled by the computer 
32. 

Fig. 3 shows a 1f common path interferometer 40. A laser 41 emits a light beam which is 
expanded by a beam expander 42 into a plane light wave of uniform intensity and directs it 

15 towards an encoder 43, i.e. a combination of an encrypted mask 43a and a complex spatial 
electromagnetic radiation modulator 43b. The light beam is transmitted through the 
encoder 43 and an image forming lens 44. The filter 45 phase shifts by 9 and optionally 
attenuates (by a factor B) the zero order diffraction part of the light modulated by the 
encoder 43. Optionally, the remaining diffraction part of the light modulated by the encoder 

20 43 may be attenuated by a factor A. The reconstructed intensity pattern is generated in the 
image plane 46 of the lens 44 and a dynamic focusing system 47 images the reconstructed 
intensity pattern Kx'.y') onto a focusing plane 48. As described for the system shown in Fig. 
1, the image l(x\y') may be detected with a camera and transmitted to a computer 49 for 
processing and authentication and the system 40 may also be controlled by the computer 

25 49. 

Fig. 4 shows details of (Fig. 4A) an off-axis read-out of an encoder 50 comprising a 
combination of a reflective mask 50a and a complex spatial electromagnetic radiation 
modulator 50b and of (Fig. 4B) an on-axis read-out of an encoder 51 comprising a 
30 combination of a reflective mask 51a and a complex spatial electromagnetic radiation 
modulator 51b with a beam splitter 52. Both configurations (Fig. 4A, Fig. 4B) may be 
utilised in the systems shown in Figs. 1-3. 
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PHASE ENCODING FOR DC PHASE FILTERING 

In the following an example of encoding an encoder and a spatial phase filter will be given 
based on a system filtering in the DC-frequency range. The exemplified system is based on 
5 a 4-f lens configuration as shown in Fig. 1 and illuminated by electromagnetic radiation in 
the visible frequency domain, hereafter simply denoted as light radiation. 

Assuming that the illuminating light is monochromatic and has a substantially flat amplitude 
profile we obtain the following spatial amplitude distribution emitted from the encoder: 



10 



15 



a(x, y) = rect — , — exp(i<Kx, y)) (1) 



yAx Ay J 



where <x(x, y) = exp(i<J>(x, y)) represent the spatially encoded phasor values and 
AxAy is the area of the input phase modulating spatial light modulator. 

It turns out to be convenient to separate a(x, y) into two terms describing a spatially 

invariant DC-value, a , and a spatially varying AC-contribution Aa(x, y). The DC-value 
can be found as: 



20 a = f f exp(i<t>(x, y))dxdy (2) 



Subsequently the AC-term is expressed by: 



Aa(x, y) = exp(i<()(x / y)) - — ^— [J exp(i<J>(x, y))dxdy (3) 



25 



The separation of a(x, y) into a spatially invariant DC-term and a spatially varying AC- 
term is an important point and will be used throughout the remaining part of this example, 
especially in the description of the spatial filtering procedure. 
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The spatial filter utilised in this example is chosen as a circular phase contrast filter 
(different transverse shapes can also be used) centred around origo in the spatial 

frequency domain, denoted by co-ordinates {f x , f y ): 



T(f r ) = 1 + (exp(i9) - l)ci 



(4) 



where f r = ^f* + f y 2 denotes radial spatial frequency and Af r describes the size of 
the circular (circ) phase filter. 



10 In the spatial frequency domain (the filtering plane) the Fourier transformation ( 3) of the 
spatially modulated light radiation from the encoder is present. The filtering operation on 
the Fourier transformed light radiation performed by the spatial phase contrast filter can be 
expressed as a simple point-by-point multiplication procedure. Subsequently the spatially 

filtered light is inverse Fourier transformed (3* 1 ) by the second Fourier lens (Fourier 
15 transformation and reflected output co-ordinates) and the resulting spatial amplitude 
distribution in the image plane (with co-ordinates (x', y') ) can accordingly be written as: 



o(x', y') = a(x', y') + (exp(i9) - l)3" 1 



3(a(x, y))circ 



f f \\ 



\ Af rJJ 



= a(x', y f ) + a(exp(i9) - l)rect — , — 

\Ax Ay 

= [exp(i<Kx', y 9 )) + a(exp(ie) - i)]rectf^- f 
L J ^ Ax Ay J 



(5) 



(X* y' 
, — 
Ax Ay; 



, one obtains: 



\c(x',yf S l + sir0||a| sin[|] - si^- - y>) + | 



(6) 
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Requiring that \o(x' 0 , y^)| 2 = o corresponding to complete darkness as the lowest 
intensity level in regions (x' 0 , y' Q ) e W Q implies: 



l + 

5 



4*1 sir (|)[H sir {|) ■ sii {*= " *° + 2 



= 0 (7) 



where the abbreviation (j> 0 = <t>(x^, y' Q ) has been used. 
The solutions to Eq. (7) are given by: 



sinf <t>- - <b + - ) ± Jsin 2 !*- - <b 0 + - ] - 1 
10 H - " 2 > \ fl )' 2 > (8, 



The requirement o < |aj < l implies that: 



sin 2 ^- - <j> e + |j = 1 => 

0 = 71 - 2(<|»- - $ 0 ) + p2n, p = 0,±1,±2, 



15 

leading to 



2 sn 



where the +sign is for 6 -values in the interval: 

20 



]i-t[ + 



(9) 



R = ^eJ => i*H< 1 (10) 



and the -sign is for 0 -values: 
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6 e - ; — + Podd 2K (12) 



The corresponding interval for (<|>- - <() 0 ) is: 

5 (»s -♦.)«]!»- J < 13 » 

Inserting the expression for JocJ , one obtains the simple intensity expression: 

\c(*, yf = 2^1 + sin(*- - <j,(x' f y>) + |J| (14) 

10 where 

JJ exp(i<|>(x, y))dxdy = AxAy|a| exp(i<|>-) (15) 

AxAy 

The phase-only transformations imply that energy is conserved: 
15 JJ |o(x\ y')| 2 dx'dy' = JJ |a(x, yfdxdy = AxAy (16) 

AxAy AxAy 

/I special case: 

The most convenient choice for a is: a = - (implying that 8 = n + p even 27i ), so that 

2 

20 the output intensity can be described as: 

\c(x' t y')\ 2 = 2[l - cos^x', y'))] (17) 

In this case the phase-*intensity mapping is described by the intervals [0; 7i] -> [0;4] . 

25 By setting a = — one obtains the following requirement to the phase function §(x, y) : 
2 
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JJ cos$(x, y))dxdy = 



AxAy 



AxAy 



JJ sin(«J»(x, y))dxdy = 0 

AxAy 



(18) 



Inserting the expression for |o(x', y')\ 2 in Eq. (16) yields: 



2 JJ [l - cos(())(x', y'))]dx'dy' = AxAy 



(19) 



AxAy 



in accordance with the first of the integral expressions in Eq. (18). 



10 



Encoding procedure: 



15 



A given intensity distribution (image) |o(x', y ')| is desired at the output side of the 
optical set-up. 

Pixellation of the image, that is generally represented in the greyscale range: 
[0; gmax] , provides the relation: 



JJ |o(x', y'fdx'dy' = AxAy => £ |o(i, j)f = ^# pi* ( 



AxAy)' 



20 • The histogram for the desired image \o(i, j)f is adjusted (adj.) within the greyscale 
range [0; gmax] , so that the previous point is fulfilled: 

\c(i tj f -Ho(i, joIU- 



• The phase values can now be calculated as: 



25 



<(>(/,./) = arccos 



gmax 
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• As before pixellation provides the relation: 
£ sin(<t>(i, j)) = 0. 

5 • The previous point can now be fulfilled by complex conjugating half the input pixels 
having the same phase value in the phase histogram. 

• The phase conjugate phase flipping provides a valuable tool (an extra degree of 
freedom) for manipulating the spatial frequency content in order to optimise the 

10 separation of low and high frequency terms at the filter plane. 

• The scheme is robust to constant phase errors across the input spatial phase 

modulator, since Eq. (14) is a function of the difference: <(>- - <t)(i, j) , only. 

Furthermore, small variations in the individual pixel phase values do not introduce any 

15 detrimental effects because the average value a , is a result of a very large phasor 
sum. 

• If the desired intensity distribution is too small to include all energy, that is, the 
histogram is scaled to maximum and the left hand side of Eq. (16) is still smaller than 

20 the right hand side, then the input phase object can be scaled until Eq. (16) is fulfilled. 
In order to obtain a scale invariant output intensity level a dynamic focusing system is 
needed. Similarly, intensity invariance can be obtained by controlling the radiated 
power from the light source. Alternatively, one can ignore the residual background 

illumination and obtain intensity levels with a gain factor of 9" (background constant 

25 equal to 1") for narrow generally shaped line structures (e.g. Eq. (6)). 

Example 1: 

A very simple example illustrating the individual steps in the above procedure will be given 
below. To simplify the example it will be considered in one dimension only. The starting 
30 point for encoding the encoder in this example is based on the following parameters: 
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a = 0.5 

9 * " (20) 
gmax = 4 



Consider the pixellated 3-step function shown in Fig. 5 to be reconstructed in the image 
plane as an intensity distribution. From the above choices of parameters one obtains the 
5 simple relation between phase values in the encoder and the image intensity values: 

|o(i)| 2 = 2[l - cosfo(i))] (21) 
To proceed from here it necessary to calculate the accumulated intensity J] |o(jl)| 2 in the 

i 

10 image to be reconstructed. The accumulated intensity is easily calculated from an image 
histogram where the x-axis represents grey level value and the y-axis represents the 

amount of pixels in the image at a given grey level value. By use of a histogram ]jT l^i)] 2 

i 

is simply found as the weighted sum of all grey level values (x-axis) multiplied by their pixel 
counting (y-axis). This describes, so to speak, the "weight" of the image. In this simple 
15 example histogram calculations are not needed since we only have 3 grey levels with well- 
defined separations. 

The value for the accumulated intensity has to obey the equality: 

20 £ H4 2 = P ix M =# P*to - 14 ( 22 ) 

i 4 

From Fig. 5 we obtain: 
J] \o(i)\ 2 = 4pixels • 0 + 4pixels • (0.5max) + 6pixels • max = 8max 



i 

25 (23) 



So that the value for max can be estimated to be: 
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max = — 
4 



(24) 



The corresponding adjusted intensity levels, |o(x)j^ d ^ , are therefore: 7/4, 7/8 and 0. These 
values can now be utilised to calculate the phase values of the encoder from the relation: 



<|>(i) = arccos 1 - 



2 > \ 



+<€» 

gmax 



- arccos 



1 - 



adj 



(25) 



where from we obtain the three phase values: 1.45 rad. 0.97 rad. and 0 rad. 
The last step needed in order to encode the encoder is that the following equality is 
10 fulfilled: 



X sin(<j>(i)) = 0 (26) 

i 

Since we have the choice to use complex conjugate phasor values (two phasors giving the 
15 same intensity level) many approaches can be taken from here. A simple approach is to flip 
every second phasor with its complex conjugate value as shown in Fig. 6. The final phase 
values used in the encoder are accordingly: ± 1.45 rad. ± 0.97 rad. and 0 rad. 



As the last step we can check whether the criteria: a = l / 2 , is actually fulfilled with the 
20 chosen phasor encoding: 

a = — (4 exp(io) + 2 exp(i0.97) + 3 exp(±L45) + 2 exp(-i0.97) + 3 exp(-£L45)) = 1/2 

14 .. 

(27) 

25 GENERAL PHASE CORRECTION PROCEDURE INTEGRATED WITH THE PHASE 
ENCODING 

In Eq. (6) we obtained an analytic relation between the phase values in the encoder and 
the resulting intensity distribution, within the region (x 1 , y') e W : 

30 
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|o(x', y')| 2 = 1 + sin(|||a| sir0) - sin(*- - y') + |j (28) 



The analysis leading to the above relation was based on the assumption that |a| is a 

constant value within the 9t ' -domain. In other words, the following approximation was 
5 applied: 



¥-1 



3(a(x, y))c ire 



\ Af rJy 



= ocrect 



Ax Ay 



(29) 



However, for certain spatial filter parameters the left-hand side of this expression will not be 
a space invariant constant value throughout the whole 9T -domain but will instead manifest 

10 slowly variations/oscillations. This will introduce small errors in the final superposition 
between the phase filtered DC-value and the direct propagated AC-signal. In order to 
circumvent this problem a technique is needed that can counteract the distortions by use of 
phase-only encoding in the components already present in the system. In what follows a 
procedure for integrating pre-distortion that counteracts the above mentioned distortions 

15 will be described that is purely based on modifying the phasor values in the encoder at the 
input side of the system. The method can also counteract other types of distortions 
inherent in a practical implementation of the system. Furthermore, the method can be 
applied in systems filtering at other spatial frequencies than DC. 



20 Procedure: 

When encoding the input phase function it is helpful to have a 'reverse' equation, 
expressing the input phase distribution as a function of an adjusted (electronic) image grey- 
level distribution, i slm , addressing the input spatial light modulator: 



25 ~ 



= i + 
gmax 



4p(x', y')| sir^Jpx', y')| sing] - sin^, y>) - y>) + |]j 



(30) 

where it has been taken into account that a(x\ y') is not considered as a constant but 
manifests a smooth oscillating behaviour within the optical image domain. The maximum 
value of l Blm is denoted gmax . 



30 
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Now, one can derive a formula for the 'grey-level correction' AI slm (x , t y') that one needs 
to apply in order to encode a phase function that compensates for the spatial variation of 
the average phase value <x(x', y') : 



4Igia(X, ' y ' ) ^ 1 + y')| "r(|)[px', y')| slag) - sin[(|)-(x', y") - y') + |)] 

r x . sfcimC*'/ yQ + AJ Jia (x , / y')) i 



gmax 
<K X '' y') = acos 



(31) 



where the second relation has been derived from the first by setting a =— and 0 = n . 

2 



By inserting the second relation in the first expression one gets: 

10 



f \ 

Aijx', y*) = 



r - 1 X slro (x<, y') - ,_ gJnaX | { |a(x', y>)\ - i] 



(32) 

This formula is however not directly useful because it is related to the histogram adjusted 
15 grey-level distribution denoted by I sla . 

One needs a formula that relates the above correction term to the 'original' input grey-level 

distribution l(x, y) that has not been modified by histogram adjustments. This is 

important since the effect of the grey-level corrections also have to be incorporated in the 
20 procedure of histogram adjustments. 

The histogram scaling gives: 

l(x, y) = -is- I sla (x, y) (33) 

slm, max 

25 

where i m and l slin max are the maximum grey-level values occurring in the original and 
the adjusted electronic grey-level distributions respectively. 
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Similarly, one can apply this relation to the intensity correction term AI 8lm and obtain: 



l(x, y) = J(x, y) + Al(x, y) = — (l slm (x, y) + Al sim (x, y)) 



Sim, max 



(34) 



5 resulting in: 



I(x, y) = -= 



2|a(x, y)| 



f if- 
I(x, y) - gmax — |a(x, y)| - - 



(35) 



In order to have 'enough dynamic range 1 in grey-levels for the correction term one can 
10 derive an inequality from the above relation by using the fact that I < gmax : 



2|(Xr 



in 



gmax I|nax f|a min | - - 

•^slm, max V J 



< gmax 



(36) 



or 



2|a 



min • gmax 



max ~~ / 



l - 



^slm, max V ' ) 



(37) 



15 



Since the first term is the dominating term in the expression for the intensity correction it 
will in practice be sufficient just to have the much simpler corrections: 



X(x,y) = -g^ 



2|a(x, y)| 

< 2|a min * CJIR3LX 



(38) 



20 



BINARY PHASE IMAGE ENCRYPTION AND DECRYPTION METHOD 



25 



Fig. 7 illustrates a binary phase image encryption method and Fig. 8 illustrates a binary 
phase image decryption method according to the present invention. In the illustrated 
example the modulating phase values are 0 or n. In Figs. 7 and 8 bright resolution 
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elements 100 have the phase shift value 0 and the dark resolution elements 101 have the 
phase shift value n. 

The encryption method comprises the steps of 

1) calculating, according to the mathematical method described above and in WO 
96/34307, the modulating phase values (0,7i) of the resolution elements forming a 
binary encoded phase pattern 102 that is imaged onto the intensity pattern 103, e.g. by 
one of the systems illustrated in Figs. 1-3, 

2) generating a uniform random distribution 104 of binary phase values (0,n), 

3) forming an encrypted phase mask 105, 4a, 23a, 43a, 50a, 51a by adding, for each 
resolution element 106 of the encrypted mask 105, the calculated phase value 107 (0 

15 or n) and the corresponding generated random phase value 108 (0 or n) of the 

resolution element, whereby a uniform random distribution of binary phase values (0,7t) 
is generated in the encrypted phase mask 105 encrypting the image 103, and 

the decryption method comprises the step of reconstructing the image 103 by 

4) adding the binary phase shift value (0,7i) of each resolution element of the encrypted 
phase mask 105 to a binary phase shift value (0,7i) of a corresponding resolution 
element of a decrypting phase mask 109, 4b, 23b, 43b, 50b, 51b with resolution 
elements having the respective calculated phase values 104, whereby the original 
encoded phase pattern 102 is reconstructed. 

The decrypting step 4) may be performed by aligning the encrypted phase mask 105 and 
the decrypting phase mask 109 in one of the optical systems illustrated in Figs. 1-3, 
whereby the original image 103 is reconstructed, e.g. for detection by a camera and 
30 recognition by a computer. 



It should be noted that in the above example, encrypting and decrypting phase values of 
corresponding resolution elements are identical. 
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Utilisation of 0 and n as the binary phase values provides a particular robust encryption 
and decryption approach with a low sensitivity to perturbations of the masks because of the 
relatively large phase difference between the binary values of the mask. 

5 In the common path interferometer illustrated in Fig. 9, the decrypting mask 120b is a 
spatial light modulator with a modulating region that is larger than the encrypted phase 
mask 120a. Thus, the encrypted phase mask 120a and the decrypting mask 120b can be 
aligned electronically by appropriate control of the spatial light modulator 120b. For 
example, the encrypted mask may contain regions with non-encrypted patterns that are 
10 imaged onto the CCD detector array 121 for the purpose of alignment. The positions of the 
imaged alignment patterns are detected and the spatial light modulator 120b is 
subsequently controlled accordingly whereby the decrypting resolution elements of the light 
modulator 120b are aligned with the resolution elements of the encrypted phase mask 
120a. 
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CLAIMS 

1. A method of decryption of an encrypted image having a non-encrypted image intensity 
pattern l(x\y') and 

5 

Encoded into a mask having a plurality of mask resolution elements (x m ,y m ) with an 
encoded phase value <|>(x mi y m ) and an encoded amplitude value a(x m ,y m ) f and 

Encrypted by addition of an encrypting phase value <J> c (x m ,y m ) to the encoded phase values 
10 <|>(x mi y m ) and by multiplication of an encrypting amplitude value a c (x m ,y m ) with the encoded 
amplitude value a(x m ,y m ), 

Each mask resolution element (x m ,y m ) modulating the phase and the amplitude of 
electromagnetic radiation incident upon it with the complex value 
15 a(x m ,y m )a c (x mj y^ and 

The method comprising the steps of 

Radiating electromagnetic radiation towards the mask, 

20 

Inserting into the path of the electromagnetic radiation a complex spatial electromagnetic 
radiation modulator comprising modulator resolution elements (x^), each modulator 
resolution element (Xd,y d ) modulating the phase and the amplitude of electromagnetic 
radiation incident upon it with a predetermined complex value a d (x dl y d )e l * d(xdyd) l the 
25 decrypting phase value M^y*) and the decrypting amplitude value a d (Xd,yd), respectively, 
of a modulator resolution element (Xd,y d ) being substantially equal to -(t> c (x m .ym) and 
ac 1 (x m ,y m ), respectively, of a corresponding mask resolution element (x m ,y m ), and 

Imaging the mask and the electromagnetic radiation modulator onto the image having the 
30 image intensity pattern l(x\y'). 



2. A method according to claim 1 , wherein the step of imaging comprises imaging with a 
common path interferometer. 
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3. A method according to claim 1 or 2, wherein the step of imaging comprises phase 
contrast imaging. 

4. A method according to claim 2 or 3, further comprising the steps of 

5 

Fourier or Fresnel transforming electromagnetic radiation modulated by the mask and the 
complex spatial electromagnetic radiation modulator, 

Filtering the Fourier or Fresnel transformed electromagnetic radiation by 

10 

In a region of spatial frequencies comprising DC in the Fourier or Fresnel plane, 

Phase shifting with a predetermined phase shift value 6 the modulated 
electromagnetic radiation in relation to the remaining part of the 
15 electromagnetic radiation, and 

Multiplying the amplitude of the modulated electromagnetic radiation with a 
constant B, and 

20 In a region of remaining spatial frequencies in the Fourier or Fresnel plane, 

multiplying the amplitude of the modulated electromagnetic radiation with a 
constant A, 

Forming the intensity pattern by Fourier or Fresnel transforming, respectively, the phase 
25 shifted Fourier or Fresnel transformed modulated electromagnetic radiation, whereby each 
resolution element (x m ,y m ) of the mask is imaged on a corresponding resolution element 
(x\y') of the image, 

The filtering parameters A, B, 6 substantially fulfilling the equation 

30 

Kx'.y 1 ) = A 2 |a(x , ,y , )e i * (x,y) + a (BA* 1 e l9 -1)| 2 
a being the average of the complex phasors a(x,y)e l * (x ' y) . 
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5. A method according to claim 4, wherein the filtering parameters A and B substantially 
fulfil that A=1 and B=1, and wherein the absorption of the mask is substantially uniform. 

6. A method according to claim 5, wherein the phase shift value G substantially fulfils the 
5 equation 



2|sin|| 

10 

7. A method according to claim 6, wherein the phase shift 9 is substantially equal to n. 

8. A method according to any of claims 3-7, further comprising the steps of 

15 Moving the DC-part of the electromagnetic radiation to a second part of the Fourier or 
Fresnel plane, and 

Phase shifting the Fourier or Fresnel transformed modulated electromagnetic radiation at 
the second part of the Fourier or Fresnel plane by 8 in relation to the remaining part of the 
20 electromagnetic radiation. 

9. A method according to claim 8, wherein the step of moving the DC-part of the 
electromagnetic radiation comprises utilisation of an optical component, such as a grating, 
a prism, etc, with an appropriate carrier frequency. 

25 

10. A method according to any of claims 3-9, further comprising the step of phase shifting 
at selected spatial frequencies constituting a region that is shaped to match the spatial 
frequency content of the phasors e j * (x,y) . 

30 1 1 . A method according to any of claims 3-1 0, wherein the step of filtering comprises 
utilisation of a spatial light modulator. 

12. A method according to any of claims 3-1 1 , further comprising the step of encoding the 
optical function of an output lens into the filter. 
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13. A method according to any of claims 3-12, wherein the step of radiating 
electromagnetic radiation comprises radiation of electromagnetic radiation of different 
wavelengths corresponding to three different colours, such as red, green and blue, for 

5 generation of intensity patterns of arbitrary colours. 

14. A method of encryption of an image having an intensity pattern l(x\y') to be decrypted 
according to any of claims 1-13, comprising the steps of 

10 Pixellating the intensity pattern l(x\y') in accordance with the disposition of resolution 
elements (x m ,y m ) of a mask, 

Encoding the mask with an encoded phase value <|>(x m ,y m ) and an encoded amplitude value 
a(x m ,y m ), and 

15 

Encrypting by addition of an encrypting phase value <t> c (x m ,y m ) to the encoded phase values 
4>(x m ,ym) and by multiplication of an encrypting amplitude value a c (x m ,y m ) with the encoded 
amplitude value a(x m ,y m ), 

20 Each mask resolution element (x mi y m ) modulating the phase and the amplitude of 
electromagnetic radiation incident upon it with the complex value 
a(x m ,y m )a c (x m ,y^ 

15. A method according to claim 14, further comprising the step of calculating the complex 
25 phasor values a(x,y)e** ( x y) of the encoder in accordance with 

l(x\y') = A 2 |a(x\y')e iwy) + a(BA* 1 e i9 -1)| 2 

For selected phase shift values 8, a being the average of the complex phasors 
30 afryje*** 

Selecting, for each resolution element, one of two phasor values which represent a 
particular grey level. 
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16. A method according to claim 15, wherein the filtering parameters A and B substantially 
fulfil that A=1 and B=1, and wherein the absorption of the mask is substantially uniform. 

17. A method according to claim 16, wherein the step of calculating the phasor values 
5 comprises 

Setting the reconstructed intensity of at least one resolution element (x 0 \ yo') of the 
intensity pattern to zero, and 

10 Calculating the phasor values e^ (x,y) of the mask in accordance with 



18. A method according to claim 17, further comprising the step of selecting the phase shift 
9 = n, selecting |a| = Vz t and calculating the phasor values e^ (x,y) of the encoder in 
accordance with 




15 



I(x\?) = 2 l + sin(<l>- -((>(*',/ 




For selected phase shift values $,<}>- being the phase of a . 



20 



25 



l(x\yO = 2(1^os(<Kx',y'))) 



J J sin(<|) (x, y)dxdy = 0 . 

encoder 
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19. A method according to any of claims 14-18, further comprising the step of encoding the 
function of an optical component, such as a grating, a prism, etc, with an appropriate 
carrier frequency, into the mask. 

5 20. A method according to any of claims 14-19, further comprising the step of adjusting the 
modulus of the Fourier transform of the complex phasors a(x,y)e'* (x,y) at specific spatial 
frequencies in order to control the range of intensity levels of the reconstructed intensity 
pattern. 

10 21 . A method according to claim 20, wherein the step of adjusting the modulus of the 
Fourier transform of the complex phasors a(x,y)e i * (xy> at specific spatial frequencies 
comprises at least one of the following measures: 

a) Adjusting the individual complex phasors a(x,y)e l<Kx y) of the resolution elements of the 
15 mask maintaining prescribed relative intensity levels between intensities of resolution 

elements of the intensity pattern, 

b) Adjusting the individual complex phasors a(x,y)e l<Hx,y) of the resolution elements of the 
mask by histogram techniques, 

20 

c) Spatially scaling the complex phasor a(x,y)e i4>(x y) pattern of the mask, and 

d) Utilising half tone coding techniques. 

25 22. A method according to any of claims 14-21, wherein each complex phasor a{x t y)e MKy) 
of the mask is selected from a set of two determined phasors with complementary complex 
phasor values a(x,y)e i * 1(x,y) and a(x,y)e i,|>2(xy) in such a way that a specific spatial frequency 
distribution of the intensity of the electromagnetic radiation in the Fourier or Fresnel plane 
is attained. 

30 

23. A method according to claim 22, wherein the phase <J>(x,y) of complex phasors 
a(x,y)e^ (x y) of adjacent resolution elements alternates between the two possible 
complementary complex phasor values a(x,y)e* 1(x * y) and afr.yje 1 * 20 ^. 
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24. A method according to claim 22 or 23, wherein the complex phasors a(x l y)e i * 1(x ' y) and 
a(x,y)e ! * 2(x,y) are complex conjugated. 

25. A method according to any of claims 14-24, further comprising the step of encoding the 
5 optical function of a Fourier-transforming lens into the complex phasors a(x,y)e iMx,y) of the 

encoder. 

26. A method according to any of claims 14-25, wherein each of the encrypting phase 
values <Mx mj y m ) is substantially equal to a value selected from a set consisting of two 

10 phase values. 

27. A method according to claims 26, wherein each of the encrypting phase values 
<j>c(x m ,ym) is substantially equal to a value selected from the set consisting of 0 and n. 

15 28. A decryption system for decrypting an encrypted image having a non-encrypted image 
intensity pattern Kx'.y') that has been 

Encoded into a mask having a plurality of mask resolution elements (x m ,y m ) with an 
encoded phase value <Kx m ,y m ) and an encoded amplitude value a(x m ,y m ) I and 

20 

Encrypted by addition of an encrypting phase value (t> c (x m ,y m ) to the encoded phase values 
<J>(x m ,ym) and b y multiplication of an encrypting amplitude value a c (x mi y m ) with the encoded 
amplitude value a(x m> y m ), 

25 Each mask resolution element (x m ,y m ) modulating the phase and the amplitude of 
electromagnetic radiation incident upon it with the complex value 
a(x m ,yja c (x m ,yje^ 

The system comprising 

30 

A source of electromagnetic radiation for emission of electromagnetic radiation for 
illumination of the mask, 

A complex spatial electromagnetic radiation modulator that is positioned in the path of the 
35 electromagnetic radiation and comprising modulator resolution elements (x^), each 



WO 00/02339 



39 



PCT/DK99/00331 



modulator resolution element (Xd,y d ) modulating the phase and the amplitude of 
electromagnetic radiation incident upon it with a predetermined complex value 
ad(Xd,yd)e i * d<xd,yd) , the decrypting phase value <j>d(x d ,yd) and the decrypting amplitude value 
a<i(Xd,yd), respectively, of a modulator resolution element (Xd.y*) being substantially equal to 
5 -<t>c(Xm,ym) and a c 1 (x m ,ym), respectively, of a corresponding mask resolution element (Xm.ym), 
and 

An imaging system for imaging the mask and the electromagnetic radiation modulator onto 
the image having the image intensity pattern l(x',y'). 

10 

29. A system according to claim 28, wherein the imaging system comprises a common 
path interferometer. 

30. A system according to claim 28 or 29, wherein the imaging system comprises a phase 
15 contrast imaging system. 

31. A system according to claim 29 or 30, further comprising 

Means for Fourier or Fresnel transforming the electromagnetic radiation modulated by the 
20 mask and the complex spatial electromagnetic radiation modulator and being positioned on 
a propagation axis of the modulated radiation, 

A spatial filter for filtering the Fourier or Fresnel transformed electromagnetic radiation by 

25 In a region of spatial frequencies comprising DC in the Fourier or Fresnel plane, 

Phase shifting with a predetermined phase shift value 0 the modulated 
electromagnetic radiation in relation to the remaining part of the 
electromagnetic radiation, and 

30 

Multiplying the amplitude of the modulated electromagnetic radiation with a 
constant B, and 
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In a region of remaining spatial frequencies in the Fourier or Fresnel plane, 
multiplying the amplitude of the modulated electromagnetic radiation with a 
constant A, 

5 Means for forming the intensity pattern by Fourier or Fresnel transforming, respectively, the 
phase shifted Fourier or Fresnel transformed modulated electromagnetic radiation, 
whereby each resolution element (x m ,y m ) of the mask is imaged on a corresponding 
resolution element (x\y') of the image, 

10 The filtering parameters A, B, 9 substantially fulfilling the equation 

l(x'y) = A 2 |a(x , ,y , )e i * (x '' y ' ) + a (BA'VM )| 2 

For selected phase shift values 8, a being the average of the complex phasors 
15 afx.yje*™. 

32. A system according to claim 31 , wherein the filtering parameters A and B substantially 
fulfil that A=1 and B=1, and for each (x,y) of the encoder: a(x,y)=1. 

20 33. A system according to claim 32, wherein the phase shift value 9 substantially fulfils the 
equation 



|oc|= 1 



2|sin— | 
1 2 1 



25 34. A system according to claim 33, wherein the phase shift 9 is substantially equal to n. 
35. A system according to any of claims 30-33, further comprising 



Means for moving the region of spatial frequencies comprising DC to a second part of the 
30 Fourier or Fresnel plane, and wherein 
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The spatial filter is adapted to phase shift the transformed modulated electromagnetic 
radiation at the second part of the Fourier or Fresnel plane by 9 in relation to the remaining 
part of the electromagnetic radiation. 

5 36. A system according to claim 35, wherein the means for moving the region of spatial 
frequencies comprising DC to a second part of the Fourier or Fresnel plane comprises an 
optical component, such as a grating, a prism, etc, with an appropriate carrier frequency. 

37. A system according to any of claims 31-36, wherein the spatial filter comprises a spatial 
10 light modulator. 

38. A system according to any of claims 31-37, wherein the spatial filter is adapted to 
perform the optical function of an output lens by appropriate encoding of the spatial filter. 

15 39. A system according to any of claims 28-38, wherein the source of electromagnetic 
radiation is adapted to radiate electromagnetic radiation of different wavelengths 
corresponding to three different colours, such as red, green and blue, for generation of 
intensity patterns of arbitrary colours. 

20 40. A system according to any of claims 31-39, further comprising a first and a second 
Fourier transforming lens, the mask being positioned in the front focal plane of the first 
lens, the spatial filter being positioned at the back focal plane of the first lens, and the 
second lens being positioned so that its front focal plane is positioned at the position of the 
back focal plane of the first lens. 

25 

41. A system according to any of claims 31-40, further comprising one Fourier transforming 
lens, the spatial filter being positioned at the back focal plane of the lens. 

42. A system according to any of claims 31-41 , further comprising one imaging lens, the 
30 spatial filter being positioned in the back focal plane of the lens. 

43. A system according to any of claims 31-42, further comprising a polarising beam 
splitter and a quarter wave plate and/or a phase filter reflecting electromagnetic radiation 
incident upon it. 

35 
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44. A system according to any of claims 31-43, wherein the spatial filter changes the phase 
of the radiation in the region of spatial frequencies comprising DC and leaves the phase of 
the remaining part of the radiation unchanged. 

5 45. A system according to any of claims 31-43, wherein the spatial filter do not change the 
phase of the radiation in the region of spatial frequencies comprising DC and changes the 
phase of the remaining part of the radiation. 

46. A system according to any of claims 31-43, wherein the spatial filter blocks the 

10 radiation at the region of spatial frequencies comprising DC and leaves the remaining part 
of the radiation unchanged. 

47. A system according to any of claims 31-46, wherein the source of electromagnetic 
radiation is a Laser. 

15 - 

48. A system according to any of claims 28-47, wherein each of the encrypting phase 
values <|>c(x mi y m ) is substantially equal to a value selected from a set consisting of two 
phase values. 

20 49. A method according to claims 48, wherein each of the encrypting phase values 
<!>c(x mi y m ) is substantially equal to a value selected from the set consisting of 0 and n. 
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Fig. 1 



WO 00/02339 



2/9 



PCT/DK99/00331 




Fig. 2 
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Fig. 4 



WO 00/02339 



PCT/DK99/00331 



5/9 




Fig. 5 
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Fig. 6 
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